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^ i ABSTRACT 

m ■ 

C\| I Is pulsar make up of strange matter? The magnetic field decay of a pulsar 

may be able to give us an answer. Since Cooper pairing of quarks occurs inside a 
sufficiently cold strange star, the strange stellar core is superconducting. In order 

^ ; to compensate the effect of rotation, different superconducting species inside a 

rotating strange star try to set up different values of London fields. Thus, we have 

O a frustrated system. Using Ginzburg-Landau formalism, I solved the problem of 

On ■ 
^ : 

Q_i' number density of vortex bundles is directly proportional to the angular speed 

O ' of the star. Since it is energetically favorable for the vortex bundles to pin to 

S-H ' 

magnetic flux tubes, the rotational dynamics and magnetic evolution of a strange 
st ar are coupled together, leading to the magnetic flux expulsion as the star slows 
down. I investigate this effect numerically and find that the characteristic field 
^ ■ decay time is much less than 20 Myr in all reasonable parameter region. On the 

c$ ■ other hand, the characteristic magnetic field decay time for pulsars is > 20 Myr. 

Thus, my finding cast doubt on the hypothesis that pulsars are strange stars. 
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rotating a superconducting strange star: Instead of setting up a global London 
field, vortex bundles carrying localized magnetic fields are formed. Moreover, the 
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1. Introduction 

What is the most stable form of baryonic matter at high density? This question is 
of both observational and theoretical interests. With the discovery of pulsars and their 
identification with the neutron stars in the late 60s, many people thought that neutron star 
matter is the most stable form of cold condensed matter at high density. This believe was 
later challenged by Witten. He proposed that at sufficiently high density, deconfinement of 
nucleons occurs. Moreover, some of the d and u quarks on of the Fermi surface is converted 
to strange quarks via weak interaction. Thus, the most stable form of matter at ultra- 
high density is a degenerate mixture of u, d, and s quarks together with a small amount of 
electrons so as to maintain overall charge neutrality flWitten 1984 ). This kind of material 



state is called strange matter. His idea was examined in detail later by Farhi and Jaffe 
(1984). 

The possibility of self gravitating strange matter, namely, a strange star, has also been 
investigated. The equation of state (EOS) calculations of strange star suggests that they 
are indeed stable objects in certain parameter range ( |Alcock et al. 1986| Haensel et al.| 



1986| ; |Benvenuto and Horvath 1989| see also Madsen and Haensel 1991 ; and Penvenuto et 



al. 19916]) . The mass and radius of a stable strange star are found to be similar to those of 
a neutron star. Consequently, some authors suggested that pulsars are in fact strange stars 
(see for example [Benvenuto et al. 1991 b| ). In principle, the question of whether pulsars are 
neutron or strange stars can be answered by comparing the core density of a neutron star 
to the strange matter transition density. Unfortunately, all the nuclear physics calculations 
to date do not yield a definitive answer. 

There are a number of difficulties in explaining some of the observational properties of 
pulsars using the strange star model. The first, and perhaps the most serious difficulty, is 
the possibility of a strange star glitch. Early EOS calculations suggested that strange stars 
have very thin nuclear matter crust. The ratio of inertial moment of the nuclear matter 
crust to that of the whole strange star is of the order of 10~ 5 . Moreover, the density of 
the nuclear matter crust is not high enough for neutron drip ( |Alcock et al. 1986| ; |Alcock 



|1991|) . In contrast, pulsar glitch observations tell us that the ratio of inertial moment of 
pinned neutron superfluid crust to that of the whole star is about 10~ 2 flAtpar 1987| ; [Link et 



|al. 1992j ; [Alpar et al. 1993|) . Nevertheless, by taking into account the existence of strange 



matter bound states, Benvenuto et al. (1991ab) showed that strange star may support a 
"strange matter crust" thick enough to account for the observed pulsar glitches. Besides, a 
more recent calculation suggested that a strange star may be able to retain a reasonable size 
of nuclear matter crust by accretion as well QBenvenuto et al. 1994j ). 

The second difficulty is the possibility of type I X-ray bursts on a strange star surface. It 
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is commonly believed that type I X-ray burst involves sudden thermo-nuclear runaway at the 
surface of an accreting neutron star (see for example |Lewin et al. 1993| ). The same bursting 
behavior cannot occur on a bare strange star surface because nuclear fuel will dissociate 
into constituent quarks immediately ( Jones 1986 ). The possibility of X-ray burst, therefore, 



requires the existence of a nuclear matter crust over a strange star. EOS calculations showed 
that the nuclear matter crust is separated from the interior strange matter by a electrostatic 
gap of a few fermi thick, thereby preventing the nuclear matter from converting into strange 
matter flAlcock et al. 1986| ; Penvenuto et al. 1991a|) . A strange pulsar magnetosphere can 



also be formed, giving rise to the observed radio pulsation ([Alcock et al. 1986| ; |Benvenuto 
et al. 1991 b| ). 



The final difficulty is the magnetic field decay time. Assuming a normal strange matter 
core, Jones (1988) pointed that magnetic field inside a strange star will not decay in Hubble's 
time. Moreover, there is no obvious mechanism for the star to retain a small residual field 
after say 10 9 years. Therefore, it is inconsistent with both the pulsar magnetic field decay 
hypothesis, and the spun- up formation scenario of milli-second pulsars QJones 19881) . Never- 



theless, the statistical evidence for pulsar magnetic field decay over a period of some 10 9 yr 
is still inconclusive. Also, milli-second pulsars may rather form by accretion induced collapse 
of white dwarfs (see for example [Bhattacharya and Van den Heuvel 1991| and [Bhattacharyaj 



|and Srinivasan 1995 for discussions on neutron star magnetic field decay). Therefore, the 



objection of Jones may not be completely well founded. Moreover, as I shall discuss in §0, 
the core of a strange star is likely to be superconducting. In this case, the magnetic field 
decay time estimated by Jones (1988) has to be modified. The question of magnetic field 
decay have been brought up in the review paper by Bailin and Love (1984). They briefly 
mentioned that the low electron density in the strange stellar core may lead to rapid flux 
expulsion. However, they did not provide a detail calculation. A major objective of this 
paper, therefore, is to perform such a calculation, taking into account the coupling between 
magnetic evolution and rotational dynamics of the star in the presence of collective effects 
such as clumping of flux tubes. 

At this moment, one has no doubt that the standard neutron star model for pulsars is 
well tested by numerous observational data. However, we still have to keep a close look at 
the alternative hypothesis, namely, the strange pulsar model ( Penvenuto et al. 1991b|) . In 
particular, we like to explore various physical properties of a neutron star and a strange star, 
and to see if there are further ways to test if which of the two models is a better candidate 
in explaining pulsar observations. 



As I shall discuss in §|2|, the core of a strange star is likely to be superconducting ( pailin 



and Love 1982| , 1984). However, this cast a problem on the rotation of a strange star. In 
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order to rotate an object with one superconducting species, a uniform magnetic field 

_ 2m*c _ , , 

b = — -n, 1 

q* 

called the London field, is set up in the object (Paym 1988|) , where m* and q* are the effective 
mass and charge of the Cooper pair, and Q is the rotational angular velocity of the object. 
But for a strange star, u, d, and s quarks are all superconducting. These superconducting 
species require different values of B to set up a rotation. Thus, we have a frustrated system. 
The situation is further complicated by the fact that all three superconducting quark flavors 
interact strongly with each other. In §g, I tackle this rotation problem using Ginzburg- 
Landau formalism. Instead of setting up a uniform London field, I find that vortices and 
localized magnetic fields are created when the superconducting quarks rotate together. Then 
in §|4], I point out that the vortices will inter-pin with the magnetic flux tubes, which alter 
the magnetic field evolution and rotation of a strange star dramatically. In particular, 
it suggests that the (superconducting) strange pulsar hypothesis is inconsistent with the 
observed magnetic field decay time in pulsars. Finally, a conclusion is drawn in §|5|. 



2. Strange Matter Superconductivity 

Using a relativistic treatment of the BCS theory, Bailin and Love (1982, 1984) suggested 
that strange matter turns superconducting at low temperatures. Using perturbative QCD 
at the one gluon exchange level, they showed that the pairing of quarks is most likely to 
occur in both the ud-du and the ss channels. The ss pairing is expected to have a gap 
matrix transforming as a color 3 and having J p = 1 + (and hence the pairing is in p- 
wave). On the other hand, pairing in the u-d system is expected to occur in the isoscalar 
channel, with a gap matrix transforming as a color 3 and having J p = + (and hence the 
pairing is in s-wavc) ([Bailin and Love 1984 ). The superconducting transition temperature is 



about 400 keV. Incidentally, the transition temperatures for neutron superfluid and proton 
superconductor in a neutron star core are of the same order of magnitude. Thus, about 
1000 yr after its supernova birth, the interior of a typical strange star is already cold enough 
for quark superconductivity ( [Benvenuto and Vucetich 1991| ). However, one should notice 
that the above estimation depends quite sensitively on the quark-gluon coupling, and may 
also change if one goes beyond the one gluon exchange calculations. Thus, the existence of 
a superconducting core in a strange star is not completely conclusive even though it is very 
likely. Nonetheless, I assume the existence of quark superconductivity in a strange star core 
in this paper. 



The quark superconductor is likely to be marginally type I with a zero temperature 
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critical field B c about 10 16 G ( |Bailin and Love 1984Q . Incidentally, the lower critical field 
of proton superconductor at the core of a neutron star is also of this order. The typical 
magnetic field of a canonical pulsar is about 10 12 G, so naively one would expect complete 
flux expulsion from the superconducting strange matter core. However, the huge electrical 
conductivity of the normal state opposes the motion of the flux, thereby, leading to a long 
flux expulsion time. In the case of a neutron star core, the flux expulsion time is ~ 10 8 years 
( [Baym et al. 1969| ). This time is shorter for strange stars because of their low electron 
density and strong quark-quark scattering amplitude. Bailin and Love (1982) predicted 
that it may be as little as 10 4 yr. However, their estimate did not take into account the 
possibility of flux clumping, which may greatly increase the flux expulsion time. Moreover, 
if the flux expulsion rate is so fast, then a Meissner state is formed in an old pulsar. As 
we shall point out in §[4.4|, the thin crustal nuclear matter is not strong enough to support 



the magnetic stress and tension of the expelled field. Therefore, we believe that during the 
life of a canonical pulsar, a meta-stable state with magnetic field penetrating through the 
superconducting strange star core has to be formed. We shall return to this point later in 
§§ 

Since the quark superfluid also couples to the color vector gluons, color superconductiv- 
ity may also be observed. Because of color confinement outside the stellar core, we do not 
have the color analogue of the ambient 10 12 G magnetic field ( Bailin and Love 1984j) . Thus, 



the observational consequences of color superconductivity remains unclear. 



3. Rotating Strange Star — Ginzburg-Landau Approach 

As I have discussed in §0J London magnetic field is set up when a single species of 
superconducting sample is rotated ( |Baym 1988| ; |heggett 199 1| ) . In fact, the presence of 



London field has been observed in terrestrial superconductors ( |Gabrera 1987| ). In § ^.1| , we 
re-derive the expression in Eq. (|]) using Ginzburg-Landau formalism as a warm up to the 
multi-superconducting-component situation like that in a strange star. In the derivation, I 
find that there is a critical angular velocity above which the uniform London field state is not 
the most energetically stable configuration. This critical angular velocity, as far as I know, 
has not been investigated in the literature. Unfortunately, such a critical angular velocity is 
too high to be attained experimentally. 
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3.1. The Case When There Is Only One Superconducting Species 



We assume that the Cooper pairs in the sample are in the s-state. Thus, the order 
parameter for the only superconducting species * is complex. In case the Cooper pairs are 



in higher angular momentum state, the order parameter will be a complex matrix flVollhardt 



and Wolfle 1990| ). We, however, shall only stick to the simple s-wave pairing derivation. 
Similar results also apply to higher angular momentum pairing states although the derivation 
becomes rather complicated. For a sufficiently small angular velocity Q and sample size, it 
is reasonable that the speed of the superfluid is not high at any point and hence a non- 
relativistic treatment will suffice. In the co-rotating frame, the Ginzburg-Landau energy 
functional reads as F = J f dV where 



/ 



-a\^\ + -1*1 + 



2m* 



h —A - m*fl xr|$ 

% c 



+ ^|VxA| 2 . (2) 



Here, m* and q* are the effective mass and charge of the Cooper pairs. Moreover, a, b > 
so that the superconducting state is prefered over the normal state. In general, a, and b are 
temperature dependent. However, a few hundred years after its supernova birth the cooling 
timescale for a strange star is much longer than its spin-down timescale. Thus, as far as 
rotational dynamics of the star is concerned, we may assume that a and b are constants. 

We consider the case when there is no external magnetic field. The presence of an 
external magnetic field may simply create extra fluxoids in the sample, which is not very 
interesting. There are two ways to minimize the kinetic energy term in Eq. (0). First, we 
can set up a magnetic field so that the ft x r contribution is cancelled by A. Alternatively, 
we can create normal cores (i.e. line defects in the form of vortices similar to that of a 
rotating superfluid). However, there are prices to pay for both cases: setting up magnetic 
field requires magnetic energy; setting up vortex needs both a kinetic energy near the vortex 
and a nucleation energy for creation of the vortex core. What really happens in the sample, 
in general, is that a uniform magnetic field together with an array of vortices will be formed 
when it is rotated. Suppose a uniform magnetic field B (parallel to ft) is set up in the star, 
then A = B x r/2. Then in order to minimize the kinetic energy term in Eq. (Q), vortices 
(or anti- vortices) have to be formed just like a rotating superfluid (see for example |Sauls 
|1989 ). Unlike fluxoids, vortices carries circulation but not magnetic field. The number of 
vortex (or anti- vortex) per unit area required is given by 



n(B) = - 



n 



2m*c 



-B 



(3) 



where k = h/m* is the circulation quantum. Obviously, vortices forms an Abrikosov lattice 
in the absence of spatial inhomogeneity. The Ginzburg-Landau free energy per unit volume 
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of the system is, therefore, given by flSauls 1989 ; Vollhardt and Wolfic 1990| ) 



F 
V 



8tt 



47T 



In 



R c 



(4) 



where p s = m^*^ is the density of superconducting species, £ is the coherence length of the 
vortex, R c is the upper cutoff length which is of the order of the inter-vortex spacing, A is 
the superconducting gap energy, and N(0) is the density of state for one spin projection on 
the Fermi surface which is given by 



N(0) 



k 



dk 



k=kp 



2n 2 h 3 c 3 



(5) 



in the extreme relativistic limit. 



Physically, the first term in Eq. (Q) is the magnetic energy, and the second term is the 
kinetic and nucleation energies of the vortices. By minimizing Eq. (f|), we find that the 
magnetic field generated in the star is 



B 



2m*c 



m ck 





n 

3 
4 



27rTW(0)A' 



if Q < Q r 



otherwise 



(6) 



where we have neglected the dependence of In R c on B because it is rather weak. The critical 
angular velocity fl c above which vortices begin to form and the translational symmetry of 
the system is spontaneously broken can be determined from Eqs. ([!]) and And it is given 
by 



2m* 3 c 2 



(7) 



Since magnetic energy scales quadratically with magnetic field while the energy required 
to form an array of vortices scales linearly with it, formation of vortices is favored at high 
rotational rates. For a typical terrestrial superconductor (high or low T c ), the effective charge 
and mass of the Cooper pairs are given by q* = — 2e and m* = 2m e respectively. By putting 
p s ~ 10~ 6 gem™ 3 , we obtain Q c ~ 6 x 10 8 rads -1 . Similarly, for neutron star matter, q* = 2e, 
m* = 2m p , and p s ~ 2 x 10 12 gem -3 , Q c ~ 2 x 10 17 rads -1 . Q c for strange star matter is 
similar to that of a neutron star. These angular velocities, unfortunately, are so high that 
the sample or the star will break apart well before the formation of vortices. In other words, 
we can neglect the magnetic energy contribution in the Ginzburg-Landau free energy all the 
time. 
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3.2. The Case When There Are Multiple Non-interacting Superconducting 

Species 

Now we examine the case when there are more than one superconducting species. We 
consider the simple case when they are non-interacting. Thus, the only coupling between 
them is that they response to the same vector potential A. If we label each superconducting 
species by j, the Ginzburg-Landau free energy density is 



2m* 



+ 



q \ 
-*-A - m*M xr $i 



— IV x A| 2 . (8) 



Once again, we expect each superconducting species to form vortices in addition to the 
generation of a common magnetic field. (Compare with Mendell and Lindblom (1991) for a 



similar study using hydrodynamics.) As we have shown in § [3.1| , we can neglect the magnetic 
energy contribution. As a result, for a uniform magnetic field B, the Ginzburg-Landau free 
energy per unit volume becomes 




*t2 



+ 



2h 



^(0)A= 



n 



2m* c 



-B 



(9) 



As shown in Fig. [I], for a fix Q, the free energy density contribution of each superconducting 
species assumes its minimum value when B equals its London magnetic field. Neglecting 
the weak dependence of lni? c on B, F/V is a continuous and piecewise linear function of B 
with "vertices" locate at the points where B equals the London magnetic field of any one of 
its superconducting species. Thus, the actual magnetic field B set up in the sample always 
equals to a London magnetic field of a particular superconducting species. In other words, 
at least one superconducting species rotates without creating vortices. While other species 
with a different London field have to create vortices (or anti-vortices) in order to rotate with 
the same angular velocity. This finding is somewhat unexpected because one would naively 
think that the value of the magnetic field set up in the sample would be a compromise 
between different London fields and all superconducting species would form vortices. As 
the superconducting species do not interact, vortices from each species will form Abrikosov 
lattices on their own (in the absence of spatial inhomogeneity) . 

What we can predict if the quarks in strange stars were non-interacting? Since the 
masses (let it be bare or consitutent) of d and u are about the same, which are both much 
less than that of s, Eq. (|l|) tells us that their corresponding London fields satisfy 

B 8 < < B ud (10) 

when Q > 0. (An exactly reverse relationship is true when Q < 0.) Since p u /m* u pa pdjm* d 3> 
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Ps/ m * s flAlcock 199 1|) , Eq. (|8]) tells us that g ud ^> g s where 



9j = for 3 = ud ' s ' 

Combining Eqs. (§), (]T0|) and (|TT[) , the slope M of our piecewise linear free energy density 
curve F/V is given by 

g s + g u d > for B > B ud 
M { g s -g ud < ioiB s <B<B ud , (12) 
~9a - 9ud < for B < B s 

where we have neglected the small nucleation energy contribution in the above estimate. 

Consequently, F/V finds its minimum when B = B ud . (Readers can verify that this 
conclusion is independent of the sign of Q.) So if different flavor quarks were non-interacting, 
only s quarks would form vortices upon rotation. 



3.3. The Case When There Are Multiple Interacting Superconducting Species 



Finally, we consider the (realistic) case when the quarks are interacting. The form of 
the free energy density must be invariant under global phase changes in any of the order 
parameters. If we consider terms up to second order in gradients and quartic in order 
parameters, the most general form of / is given by (compare with Alpar et al. 1984b where 
they have omitted the quartic self-interaction terms) 

/ = ^|VxA| 2 + W-^| 2 + ^ 

3 \ j ) j,k 

^u v ^j) ■ wtk)*m} , (13) 

where b is a real symmetric matrix, /i is an hermitian matrix, and 

Vj = — + ^A - m*n x r (14) 

i c J 

is the covariant momentum operator in the co-rotating frame for the j-th species. Since 
the free energy has to be bounded from below, we further require all the eigenvalues of the 
symmetric matrix b to be positive. 
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To explicitly show the drag effect of one superconducting species on the other, we write 
the order parameters as |^-| exp(iifj), where ipj are the phases of the order parameters. 
Then, the velocity of the superconducting species j in the rotating frame is given by 

1 h q* 

v,- = — P v ^ 7 « — Vw,' + -^-A - U x r , (15) 

m*^j J 3 m* 3 m*c 

where we have neglected the spatial variation of |^-|. From Eqs. ([Tj^) and (|T5j), we have 

/ = -L| V x A| 2 - £ + \H [b jk \^j\W + PjkVj ■ v fc ] , (16) 

where 

Pii = ™jl*il 2 + /%m* 2 |*;| 4 + 2^mf|*/ (17a) 

and 

p jk = (Pjk + 2vjk) m*m*|^f |^ fc | 2 (17b) 
for j ^ fc. Thus, once a superconducting species moves, Eq. ( |I6"D tells us that in general it is 
energetically favorable for the other superconducting species to move along with it as well. 
A similar conclusion can be reached by using three-velocity hydrodynamics ( [Khalatnikov 
n?57 ; |Andreev and Bashkin 1976] ; |Vardanyan and Sedrakyan 198 1| ) . 



At the strange matter density, which is about 5 x 10 14 gmcm 3 , the "fine structure 
constant" for strong force a s is about 0.5 - 0.6 ( Penvenuto et al. 1991b|) indicating that 



QCD is the dominant interaction between the quark Cooper pairs. The values of the off- 
diagonal terms of the matrix p, which measure the strength of the drag, can be calculated (in 
principle) from the QCD interaction Hamiltonian. Although the large value of a s prevents 
us from using perturbation theory, we believe that the effective coupling constants faj and 
Vij in Eq. ( |TB| ) are at least of order of a s . Thus, the drag force between different quark 
flavors are so strong that up to first order approximation, all three quark flavors move at 
the same velocity except possibly near the superfluid cores. The co-moving approximation 
greatly simplifies our effort to find the minimum system configuration of a strange star. 

In order to minimize the drag energy, all the superconducting species have to share 
common normal cores. We call such a configuration a vortex bundle. The circulation of the 
j-th superconducting species around a circle of radius r centered at the axis of rotation of 
the star in the co-rotating frame is given by 

= / Vj • d\ 

' V x (SI x r) • dS 



cf—Vi Pi dl+\ [B-dS- 
J m* m*c J 

hNjL q*$ v L q*nr 2 B g 

°— + ^ h - 2nr 2 Q 



rrij rrijC rrijC 
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where f2 is the angular speed of the strange star seen by an external inertial observer, Nj is 
the number of vortex quantum per bundle for species j, L is the number of vortex bundles, 
<3>„ is the magnetic flux in the core of a vortex bundle, and B g is the global uniform (London) 
magnetic field in the star. Thus, the vortex bundle density T> is given by 

V hNj + = 2m*n - for all j. (19) 



Using the same idea, we can show that the speed of superconducting species j at a small 
distance r from the core of a vortex bundle as seen in the co-rotating frame is given by 



Vj(r) 



2Trm* r 



q*B g r 



2m* c 



Thus, the co-moving requirement of super-currents at all spatial points requires that (a) 
q*B g /2m*c — Q is a constant for all j, which is possible only if the global uniform London 
field B g is zero; and (b) hNj/m* + q*§ v /m*c is a constant for all j, implying that the integers 
Nj, and the magnetic flux <$> v are chosen in such a way that 

h -^l + ^ = K (2!) 

m* m*c 

for some constant K ^ for all j. K can be interpreted as the circulation of a vortex bundle 
in this strongly interacting superconducting system. 

The necessary and sufficient conditions for the existence of solution in Eq. ([H]) are 
proven in Appendix 0. It turns out that the co- moving constraint (Eq. ([H])) is very strin- 
gent: In general, solution may not exist for a system involving more than two species of 
superconducting Cooper pairs. And in the case the solution of Eq. (|2l|) does not exist, the 
only way out is that superconductivity in some species are destroyed. Luckily, in case of the 
strange star matter, there is only two species of Cooper pairs, namely, ud-du and ss, so that 
the solution in the strongly interacting limit exists. In this limit, we expect all the quark 
flavors to form quantized vortices when the star rotates. The normal vortex core of each 
superconducting species shares a common region of space, in the form of vortex bundles. 
Magnetic field may be present in the vortex bundle cores. Moreover, stellar rotation is made 
possible by the formation of vortex bundles rather than by a uniform London magnetic field. 

In case of the strange star matter, the ud-du and ss Cooper pairs have charges g* d = e/3 
and q* = — 2e/3, respectively. Also, the effective mass of a strange quark at the strange star 
interior m s ~ 175 MeV ( |Benvenuto et al. 1991b]) . Therefore the effective mass of an ss 



Cooper pair, m*, is approximately 350 MeV. The effective mass of a u or d quark, m u rf, is of 
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order of 10 MeV, giving m* ud ~ 20 MeV. Since (m ud , q ud ) and (m s , q s ) are linear independent, 
the solution of Eq. (^) is given by (see Appendix |A|) 



K 



- q* s N ud 



(22) 



where N ud and iV s are the number of quanta per vortex bundle for ud-du and ss Cooper 
pairs respectively, and K ^ 0. In the zero temperature limit, the system will choose the 
ground state configuration out of the above infinitely many solutions. Similar to Eq. (f|), the 
Ginzburg-Landau free energy per unit volume in the strongly interacting limit is given by 



F 
V 



E 



h 2 Psj Nf 
Amm* 2 



In 



,17. 



^•(0)A| 



+ 



8tt 2 A 2 



W\ 



(23) 



where A is the penetration depth, and T> is the number density of vortex bundles. From 
Eq. (|T§D, V(hNj + q*§ v /c) = 2m*Q in the strongly interacting limit. Consequently, 



F 
V 



E 



hNjC + q*$ v 



h 2 p SJ N. 



2-Km* 



In 



R, 



+ nm*^N j (0)A 2 j 



+ 



<& 2 v ml 



c\Q\ 



^ 2 \ 2 (hN lC + qt§, L 
increases (approximately) linearly with the angular speed \Q\ of the star. 



(24) 



Now, we estimate the average free energy density for strange star matter. BCS theory 
tells us that A = 1.76kT c for s-wave paired superconductor. Putting kT c ~ 400 keV for 
strange stellar matter ( [Bailin and Love 1984[ ), we find A ud ~ 700 keV. The number densities 
for u and d quarks are given by ( jAlcock et al. 198"S| ) 

2a s 



n ; 



71 



for j = u,d, 



(25) 



where pj is the chemical potential. Putting a s = 0.5 and \i u ~ p d ~ 400 MeV ( |Benvemito 
et al. 1991b|) , n u w n d ~ 5.7 x 10 38 cm -3 . Since the electron number density is much less 
than that of the quarks ( [Bailin and Love 1984| ; |Alcock et al. 1986| ; [Benvenuto et al. 1991b| ), 
the charge neutrality condition reads 



-n„ n d n H 

3 3 3 



n e _ 







(26) 



Therefore, we expect n s 



n d . 



For both u and d quarks, Ep pa fi ~ 400 MeV ( |Bailin and Love 1984| ), and hence from 
Eq. ([5], N ud (0) ~ 6.6 x 10 41 erg -1 cm -3 . Also, m u ~ m d ~ 10 MeV. The coherence length £j- 
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can be estimated in the framework of relativistic Ginzburg-Landau theory. Bailin and Love 
(1984) argue that 

g - 9h l a3) , (27) 

where ((3) ~ 1.20, and pf^ is the Fermi momentum of species j which is approximately equal 
to Ep./c inside a strange star. Thus, ^ ~ 8.0 fm. Since the quark superconductor is likely 
to be marginally type I, A ~ v^2£«d ~ H-3 fm. Assuming \n(R C j/^j) to be of order of 10, I 
numerically compute the average free energy density in Eq. (|24T ) for all possible values of N ud 
and N s . I find that the ground state configuration for strange star matter is achieved when 
N s = 1, N U( i = 0. Consequently, the values of K and in the ground state configuration 
are given by 

K = V„* ~ °- 012 cm2 s_1 ( 28a ) 

and 



m *s - 2m *ud 



= . 3kcm * ud . ~ 8.0 x lO" 8 G cm 2 (28b) 
e(m* - 2m* ud ) 



respectively. 

In summary, the star behaves quite differently in the non-interacting and strongly inter- 
acting limits. In the non-interacting limit, a global magnetic field is set up, and vortices are 
formed in all but one superconducting species when the system rotates. In contrast, vortex 
bundles in the form of Abrikosov lattice is present in a rotating strongly interacting system. 
Moreover, the global uniform London magnetic field is no longer present. Regarding the 
coupling between the quarks as a tuning parameter, it is interesting to map out the phase 
diagram of this system. And this will be carried out in future works. 



4. Vortex Inter-pinning And Its Astrophysical Consequences 

The formation of vortices in a rotating superconducting strange star has at least three 
possible observational consequences if one assumes that the observed pulsars are in fact 
strange stars instead of neutron stars. First, the heat capacity of the star and the cooling 
processes are modified and its effect on the strange star cooling has been studied ( Benvemlto 



and Vucetich 1991| ; [Page 199"2[ ). In this section, we concentrate on the second effect, namely, 



the magnetic field decay due to inter-pinning of vortex bundles and the magnetic fluxoids. I 
shall discuss briefly possible magnetic field alignment due to inter-pinning in §[4.2| as well. 



Nucleation energy is required in the creation of normal cores in both vortex bundles 



-14- 



and magnetic fluxoids in a strange star. Therefore, it is energetically favorable for a vortex 
bundle to "pin" to a magnetic ffuxoid so that a lesser volume of normal strange matter has 
to be nucleated. Thus, the dynamics of vortex bundles and magnetic flux tubes are coupled 
together. Similar to the case of superfluid neutron vortices in a rotating neutron star ([Alpar 



|et al. 1993p , Eq. (^4|) tells us that in the strongly interacting limit, the quark vortex bundle 
density is directly proportional to the angular speed of the rotating strange star. As a 
strange star spins down, its vortex bundles in the star core have to move radially outward 
from the rotational axis (and eventually some of them will annihilate near the stellar surface). 
The pinning of quark vortex bundles and magnetic flux tubes implies that the rotational and 
magnetic evolution of a strange star are coupled. 

Alternatively, the vortex bundle and the flux tubes can pin together by interaction 
of their core magnetic fields. When two magnetic carrying wires are placed together, an 
energy change of Bi ■ J$2V/4:1T is expected, where V is their interacting volume, and Bj are 
their magnetic field strength. Depending on their relative field orientation, the two wires 
experience either mutual attraction or repulsion. In the former case, the wires will pin to 
each other; and in the latter case, the wires will avoid each other, and hence effectively "pin" 
to the inter- wire spaces. In fact, similar ideas have been applied to the proton fluxoids and 
magnetic vortices in neutron star cores; and their possible observational consequences have 
been studied flSauls 1989| ; |Srinivasan et al. 1990| ; |Chau et al. 1992| ; |Ding et al. 1993| ) . 



4.1. Estimation Of Pinning Energy Per Intersection 

Now we estimate the pinning energy per intersection of a vortex bundle with a flux tube 
by both pinning mechanisms. The pinning energy due to nucleation per intersection is given 

by 

(^)„ U ci^^E^(0)A?^- (29) 

Z 3 

The intersection volume for species j, Vj, is given by 

Vj « avrACg , (30) 

where £j is the coherence length for species j, A/fl ux is the number of flux quantum in a flux 
tube, and a is a geometrical factor depending on the angle between the magnetic flux tube 
and the vortex, 6, and their elastic moduli. For stiff magnetic flux tubes and vortex bundles, 
a ~ 2 cosec 101. 



We estimate the nucleation pinning energy as follows: for an s-wave superconductor, 
BCS theory tells us that A = 1.76kT c ; and for a p-wave superconductor, A = 2AkT c ( |Baym| 
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and Pethick 1975). Putting kT c ~ 400 keV for strange stellar matter ( [Bailin and Love 1984Q , 
we obtain A ud ~ 700 keV and A s ~ 960 keV. For ud-du quark Cooper pairs, E F rs \i ~ 
400 MeV ( [Bailin and Love 1984[ ). And beta equilibrium requires E Fs « /x s = fi d « 400 MeV. 
The coherence lengths are given by Eq. ([?7]), giving us C, u d ~ £ s ~ 8 fm. So from Eqs. (Q), ( [29] ) 
and (|30|) , we obtain 

(E p ) nud ~ 2.4< x 2 cosec |0| MeV. (31) 
Similarly, the magnetic pinning energy is given by ( |Chau et al. 1992| Ping et al. 1993| ) 



(P\ ~ [Byortex • B fluxoid | i^ e <l>^AA/' fll { x , , , , 

(^)mag ^ ^ ^ l COt ^l » ( 32 ) 

where -B c is the critical magnetic field. Since the flux quantum for a strange star flux tube 
equals hc/2(2e/3) = 3foc/4e, and A ~ 11.3 fm, we obtain B c ~ 7.7 x 10 16 G. By putting the 
value of obtained in Eq. ( |28b|) , we find 



(£pLa g ~ eso-^fluxl cot d \ MeV - ( 33 ) 

So, unless the magnetic flux tubes and the vortices are perpendicular to each other (this 
happens when the rotational and magnetic axes of the star are orthogonal to each other), the 
magnetic pinning energy is approximately two order of magnitude larger than the nucleation 
pinning energy. Therefore, we shall use the magnetic pinning energy in our subsequent 
calculations. 



4.2. Forces Acting On Flux Tubes And Vortices 

Following Ding et al. (1993), the major forces acting on flux tubes and vortices of a 
spinning down strange star are summarized below: 

(a) Inter-pinning force: Since it is energetically favorable for the magnetic flux tubes and 
the vortex bundles to pin with each other, a force of order of E v j £ is experienced by both 
parties when one tries to pull them apart. (Similar assertion is true when the magnetic 
pinning force is repulsive.) This inter-pinning force encourages the magnetic flux tubes 
and vortex bundles to move together. 

(b) Thermal activation: Since E p is of order of 100 MeV and the interior temperature of 
the star is < 10 keV, random thermal noise can be an important source of "force". In 
particular, thermal activation can depin a vortex with a flux line. 
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Magnus force: The magnetic flux tubes couple to the thin crust of the strange star via 
electromagnetic interaction between its strong core magnetic field and the background 
electron plasma. Therefore, the flux tubes co-rotate, and hence spin down, with the 
crust of the star. The vortex bundle core, due to inter-pinning, may spin down with the 
crust as well. Eq. (^) tells us that the only way to slow down the angular speed of a 
rotating superconducting quark fluid is by reducing its vortex bundle density. This can 
in turn be achieved only by moving the vortex bundles radially outward from the axis 
of rotation of the star. So as the star slows down, a steady angular velocity difference 
between the superconducting quark fluid and the crust is developed. Consequently, a 
hydrodynamic force, called the Magnus force, acting on the vortex bundles is developed. 
Magnus force tries to push the vortex bundle radially outward from the rotational axis of 
the star, thereby reducing the angular velocity difference between the superconducting 
quark fluid and the crust. The combined effect of the Magnus force and the inter-pinning 
force may cause a spin down induced magnetic field expulsion in a strange star. The 
Magnus force also defines a preferred direction for the thermal activation (i.e., thermally 
assisted creeping). 

For a inter-winded network of pinned magnetic flux tubes and quark vortex bundles, the 
average force per unit length acting on a magnetic flux tube due to the Magnus force 
acting on the vortex bundles is given by (compare with ping et al. 1993| ) 



f/,Ma g ~ ^pKrue r , (34) 

where r is the distance from the rotational axis, u is the angular velocity difference 
between the superconductor (Q s ) and the stellar crust (f2 c ), K is the circulation given 
by Eq. ( |28aj ), p is the total matter density, Lf and L v are the total number of flux tubes 
and vortex bundles respectively in the star. In fact, Lf and L v can be determined by 
considering the total magnetic flux and the circulation of superconducting current of the 
star. They are given by 

Lf = 47re 5-f^' , (35) 



and 

^RLM _ 47r^ tar (m s - 2m ud ) 



fi| , (36) 



K h 

where -R s tar and -B st ar are the radius and total magnetic field strength of the star respec- 
tively. 

The angular velocity difference, u, is likely to remain at its steady state value, which 
can be deduced from the vortex creep theory ( [Alpar et al. 1984a| ; [Alpar et al. 1993|) and 
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is given by 



kT 

= ±^77 — rsinh" 
pK b p r\ 



lol f'-^jA 
| s 2 1 exp ' 



XpKb p r 



|cos#| {elBst 



npKr V 3/ic 



1/2 

= ±^cr , (37) 



where we have assumed that E p 3> fcT. Here t>o ~ 10 13 cms -1 is the microscopic creeping 
speed, bp ~ (1/| sin#|) x (7ri? 2 tar / ' Lf) 1 / 2 is the mean distance between successive pinning 
sites along a vortex bundle, and u) cr is the critical angular velocity lag above which it is 
no longer energetically favorable for the vortices and fluxoids to pin together. The plus 
(or minus) sign is taken if the outward moving speed of vortex bundles is greater (or 
less) than that of magnetic flux tubes. Therefore, the maximum possible average force 
per unit length vortex bundles can exert on a flux tube equals pKR staT u cr L v / Lf. 

Finally, if the spin down of the star is mainly due to its dipole radiation lost, then the 
outward moving speed of the vortex bundles is given by (ping et al. 1993"D 



S 2 t _ T .i? s+ „_f2 2 sin 2 9 ,„_,,. 
v v (t) w stai ^ , (38) 

where / is the moment of inertia of the star. 

(d) Buoyancy force: The presence of magnetic stress in the core of a flux line decreases the 
internal density of quark matter ( |Muslimov and Tsygan 1985| ; |Harvey et al. 19861 ; |Jones 



1987| ). Thus, the buoyancy force per unit length experienced by a flux line is 



l-buoy 



128e 2 7r 2 A 2 c2 r 128e 2 vr 2 A 2 iL 



(39) 



where g is the local acceleration due to gravity, c 2 = dP/dp ps gR s t&r is the squared 
sound speed. 

Tension: The combined pushing of all the vortex bundles pinned onto a given magnetic 
flux line may globally bend the flux line. In effect, a tension force is developed which tries 
to resist further deformation. The average tension per unit length is given by ( Harvey 
et al. 1986Q 



ftens = - 111 U J T e Gr w ^~ ln UJ Gr ' (40) 

where s c is the radius of curvature of the flux tube. 
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(f) Electron drag force: As a flux tube drifts out of the core, it will experience a drag force 
arising from the scattering of the degenerate relativistic electron. This drag force limits 
the speed of the flux tubes and hence the rate of magnetic field decay. In the absence 
of clumping of flux tubes, the drag force per unit length is given byQ fllVluslimov and 
Tsygan 19851 ; IHarvey et al. 198ft Uones 19871 ) 



27tt n e h 2 cNl£v } 
1024 E f (e)\ 



fdrag i no a m /_\ \ ^ r ' (^1) 



where Ef(e) is the electron Fermi energy, Vf is the velocity of the flux tube, and n e is 
the electron number density and is given by ([Alcock et al. 19861) 

- = (42 > 

The electron chemical potential in a strange star is about 20 MeV ( [Benvenuto et aE 



1991b ), giving us n e ~ 3.5 x 10 34 cm -3 . Thus, the electron number density is some two 
orders of magnitude smaller than that of a neutron star. Consequently, the electron drag 
force in a strange star is much weaker than that in a neutron star. 

The average radial velocity of a flux tube is approximately equal to its steady state 
radial velocity Vf, which can be calculated from the force balance equation 

ff,Mag + fbuoy ftens fdrag 

= . (43) 
After finding Vf, the magnetic field -B sta r of the star can be computed by solving the equation 



(compare with the f2 equation in [Alpar et al. 1984a| ) 



4 t ar = • (44) 

-'tstar 



4.3. The Alignment Torque 

In general, magnetic flux tubes incline at a non-zero angle 6 to the rotational axis of the 
star. Thus, the magnitude (and sometimes even the direction) of force acting on a flux tube 
by the vortex bundles changes as we go along the flux tube itself. The combined forces acting 
on a magnetic flux tube due to vortex bundles, therefore, produce a total force together with 



2 Note that there is an extra factor of 9/4 because the flux quantum for proton superconductor and strange 
matter superconductor are different. In addition, the case of flux clumping will be addressed in £4.5. 
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a torque (Ruderman 1991abc). While the total force may lead to a magnetic field decay in a 
spinning down strange star, the torque tends to push the field lines towards the equator as 
the star spins down. Similarly, it tends to push the field lines towards the poles as the star 
spins up ( |Rudcrman 1991 b| ). 

The effect of this torque on the time evolution of magnetic field is not completely clear. 
If the crustal material of a strange star, which makes up of nuclear instead of strange matter, 
is strong enough to support a shear, then the direction of magnetic field will not change with 
time. On the other hand, if the crustal nuclear material is weak and brittle, direction of 
magnetic field may change via a series of "crust cracking" (Ruderman 1991abc). A more 
complete discussion on the effect of this alignment torque will be reported in future works. 

4.4. The Crustal Magnetic Field 

So far, I am concentrating on the magnetic field in the core of the star. In this subsection, 
I show that the presence of crustal magnetic field does not seriously affect the magnetic field 
evolution of a strange star. The density of the strange stellar crust must be below the 
neutron drip density (~ 4.3 x 10 11 grncm" 3 ), otherwise the dripped neutrons can convert 
into strange matter by making contact with them. Since the most favorable nucleus just 
below neutron drip is 118 Kr flBaym and Pethick 1975Q , spacing between the nuclear lattice 
just below the neutron drip b is about 70 fm. Hence, the Young's modulus of this lattice 
is about (Ze) 2 /fe 4 ~ 10 29 dyncm -2 . (The Young's modulus of the nuclear lattice is much 
smaller if it is not prefect or if the density of the crust is lower than for the neutron drip.) If 
the magnetic flux is completely expelled from the core, the inner layer of the normal matter 
crust has to sustain a stress of ~ B 2 AR/47r per unit length, where AR ~ 100 m is the 
thickness of the crust. The magnetic stress lengthens the nuclear lattice; however, such a 
lengthening must not be greater than the < 10 fm electro-static gap between strange and 
nuclear matter in the star. Otherwise, strange matter conversion will take place and the 
crust will be destroyed. After some computation, I find that the maximum field and flux 
the nuclear crust can sustain are ~ 10 7 G and ~7x 10 17 Gem 2 respectively. This is much 
smaller than the initial magnetic flux of a strange pulsar. So we have two possibilities: (a) 
there is an efficient field decay mechanism operating in the thin nuclear matter crust; or (2) 
the crust breaks at some point when majority of the field in the core is expelled, eventually 
making the star a bare strange matter object. In either case, the contribution of the magnetic 
field in the crust does not play an important role in the magnetic evolution of the entire 
star. Thus, we shall neglect the presence of crustal field in our subsequent analysis. 
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4.5. Magnetic Field Decay — Numerical Results 

I numerically compute the magnetic field decay due to inter-pinning of magnetic flux 
tubes and vortex bundles. I assume that the initial magnetic field -B star (0) is uniformly 
distributed in the stellar core, and we also neglect the existence of a thin nuclear matter 
crust. The initial angular speed of the star is set to be Q(0). The effect of the alignment 
torque is also ignored. To simplify calculations, we follow Ding et al. (1993) to combine the 
buoyancy force (Eq. (j39f)) and the tension force terms (Eq. (|40|) ) because of their similar 
dependence on parameters. We write 

/buoy + /tens = 7/buoy , (45) 

where 7 = 1 — (-R star /2s c ) ln(A/£) is a time dependent quantity of order of unity in most 
parts of the star ( [Harvey et al. 1986| ). 

What happens if the magnetic flux lines form a clump? Provided that the size of 
such a clump is smaller than the electron mean free path, electrons may collectively scatter 
with a number of flux lines within a clump leading to a dramatic change in the drag force 
( [Ruderman 1992| ). In addition, clumping leads to an increase in the local density of pinning 
centers. Moreover, the force between flux tubes within a clump may be important. Thus, our 
mean field estimate of the force acting on the flux tubes by the vortex bundles will change 
as well. 

We model the clumping effect by a renormalization scheme. By coarse graining to the 
level of a clump, the behavior of the flux tubes inside a clump is similar to that of a single 
magnetic flux tube carrying the same amount of flux as the clump provided that the size of 
the clump is much smaller than the electron mean free path. In addition, the magnetic field 
strength in this single magnetic flux tube equals to that in each of the individual flux tubes 
in the clump. That is, we may replace A/fl ux by A/fl ux A/' c i um p. As a result, the electron drag 
force, vortex acting force, and the buoyancy force on a clump are given by 

/drag,coll ~ A/" clump /drag , (46a) 

/Magxoll ~ A/"cl ump // i Mag , (46b) 

and 

/buoy, coll ~ A/clump^/buoy (46c) 

respectively. So the electron drag force becomes dominant when clumping of flux lines is 
serious or when the quark superconductivity is extremely type I. And from now on, the 
symbol A/fl ux should be interpreted as the product of the number of flux quantum in each 
flux tube and the number of flux tubes in a clump. 
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Let us consider a 1.4M strange star with (core) radius 11 km whose density in the outer 
region of the strange core is around 5 x 10 14 gmcm -3 . The moment of inertia of the star is 
about 1.6 x 10 45 gmcm 2 QBenvenuto et al. 1991 bp . We also fix 9 = 45°. We take a typical 



star with initial magnetic field _B star (0) = 10 12 G, initial angular speed Q(0) = 2000 rads -1 
(and hence an initial period of about 3 ms), 7 = 0.5, and A/fl ux = 1 as our "reference" star. 
Then the effects of the initial magnetic field, initial angular speed, the value of 7 in the 
effective buoyancy force, and the number of flux quanta in a flux clump AWx on the stellar 
magnetic evolution can be studied by varying these parameters one at a time. 

We first investigate the case when there is no clumping, and the system is almost 
type II. So we set Afl ux = 1- As shown in Figs, f^a and |3|a, both the magnetic field and 
angular speed remain approximately constant for a while which are followed by power law 
decays with exponents equal —1/4. Their transition times decrease with increasing initial 
field. In addition, the radial velocities of both the vortices and the fluxoids moves as a whole 
most of the time (see Fig. ^a). 

We can explain both the field and angular speed evolution in a simple way. Since 
the initial magnetic field is not too high, the star only experiences a modest spin down. 
Thus, the vortex bundles are effectively pinned to the flux tubes (instead of thermally creep 
through them) leading to what Ding et al. (1993) called a "co-moving phase". In this phase, 
-^star(t) ~ Q{t) and hence the dipole spin down of the star is given by 

• ^ s 6 t Asin 2 g ^ i? 2 ar (0R 6 ta^ 5 sin 2 g 

3ic 3 ~ 2/c 3 fi 2 (0) ' 1 ' 

Upon integration, it is easy to show that 

'ABl^)Rt t3I tsm 2 9 , ] 



3/c 3 fi 2 (0) fi 4 (0) 



-1/4 


' t 1 


-1/4 




.to + fi 4 (0). 





(48) 



Thus, Q (and hence B Btax ) is almost a constant when t -C to and Q decays as a power law 
with an exponent —1/4 when t ^> t . This is consistent with the field and angular speed 
evolution we have plotted in Figs. 0a and ^|a. 

Now, we go on to consider the case when clumping is important, and when the quark 
superconductor is almost type I. We illustrate the situation by putting A/flux = 10 6 . As 
shown in Figs. 0b and |^b, when the initial field is low (< 3 x 10 12 G), the magnetic and 
spin evolution behave almost in same way as in the case when Afl ux = 1- The huge value of 
A/fl ux implies that the electron drag force may be dominant. In order to push the flux tubes, 
the vortex bundles need to acquire a strong Magnus force by increasing the value of the 
angular velocity lag u. Nevertheless, |o;| cannot exceed its critical value u cr . So, when the 
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star spins down too quickly, the vortex bundles in the core have no choice but to thermally 
creep through the flux tubes leading to what Ding et al. (1993) called the "forward creeping 
phase". This picture is consistent with the numerical finding that forward creeping phase 
lengthens with increasing initial magnetic field (see Figs. [|b, c and d). The small radial 
velocity of the flux tubes at this phase implies that -B s tar remains almost constant. It is until 
the onset of the co-moving phase due to a much slower spin down rate at later times that 
I? s tar and Q decay like t~ 1//4 (see Figs. ||b and||b). 

Now we discuss a more interesting case when -B star (0) = 10 13 G and AVx = 10 6 . Fig. |]b 
shows that -B s tar decays exponentially with a characteristic time of about 5 Myr. However, 
at about 100 Myr, the decay stays almost constant for a while, and then it decays further 
like t -1 / 3 . Besides, Fig. [5]b shows that f2 decreases like a power law with an exponent —1/2 
for almost 10 Myr. After that, Q stays almost constant for a while which is then followed 
by a power law decay with an exponent about —1/6. 

We can explain the behavior as follows: The rapid spin down due to high value of -B s tar 
leads to the forward creeping phase (see Fig. f|d). So, 

^ ^ Q 3 i4A r (0)sin 2 g 



and hence 



2 a i -1/2 

+ 



(50) 



31 c 3 n 2 (o)_ 

This accounts for the exponential decay in Q starting from 3/c 3 /25g tar (0)i?g tar f2 2 (0) sin 2 6 « 
5 yr. About 10 3 yr or so, the star rotates so slowly that vortex bundle density becomes 
very low. At this moment, buoyancy becomes the dominant driving force for field decay and 
hence 

• 2B stax v f _ lQ>yE f (e)cB st&r _ B star 

-Rstar 37T 3 e 2 n e A-R 2 tar A/' flux r buoy 

That is, -B s tar decays exponentially with a characteristic time Tb uoy of about 4 Myr. Finally 
at time > 10 Myr, dipole spin down of the star becomes very ineffective due to the small 
values of Q and -B star . At this time the radial velocity of flux tubes is faster than that 
of the vortex bundles and the star enters the "reverse creeping" phase (see Fig. |]d and 
Ding et al. 1993]) . The pinning force prevents the flux tubes from moving too fast and 



hence the rate of field decay is decreased. The force balance equation for the flux tubes 
reads /buo y ,coii - /M ag ,coii = /drag.coii ~ (compare with Eq. (|3|)). Since / b uo y ,coii is time 
independent, /jviag.coii oc \Qu)/B Btax \, and 10 ~ —^ cr oc l-Bstarl 1 ^ 2 , the force balance equation 
implies that -B sta r ~ As a result, Q(t) ~ t~ l / & and -B sta r(^) ~ t" 1 / 3 at the very late stage 
of the evolution of the star (see Figs. |2|b and |||b). 



-23- 



We now turn to the study of the effect of f2(0) to the magnetic evolution of the star. 
Fig. [| shows that the field decays more rapidly (due to the faster onset of the power law 
decay) as Q(0) increases when the star is in co- moving phase. This agrees with the crossover 
time estimate given in Eq. (f£8|). However, if the star is in the forward creeping phase, the 
push received by the flux tubes due to pinned vortices saturates (and thus, independent of 
Q). So in this case, the value of Q(0) has little effect on the magnetic evolution of the star. 
I have verified this in some of our runs. 

Next, we consider the effect of 7. Again, the only situation where the value of 7 can 
seriously affect the magnetic evolution of a star is when the pinning force is weak at some 
moment as compared to the buoyancy force; and this happens when both 5 star (0) and A/fl ux 
are large. Fig. ^| show the magnetic field evolution when -B star (0) = 10 13 G and A/Vx = 10 6 . 
We see that the evolution profile depends quite sensitively on 7. When 7 = 0, we observe 
the -B s tar ~ i -1 / 4 behavior at late times, indicating that the star is in the co-moving phase. 
This is reasonable since the only force that push the flux tubes out from the core comes 
from the vortex bundles. But when 7 > 0.5, buoyancy force becomes the dominant. Thus, 
the magnetic evolution follows an exponential decay when t < 10 6 yr and -B s tar ~ t -1 / 3 at 
t ^ 10 9 yr, indicating that the presence of forward creeping and reverse creeping phase at 
early and late times respectively. 

Finally, we consider the effect of A/fl ux - For B stai (0) = 10 12 G, Fig. [7|a tells us that the 
behavior of the star does not have any visible change A/fl U x increases from 1 to 10 3 . They are 
all in the co-moving phase. This observation is consistent with the fact that to i n Eq. fl48|) 
is independent of A/fl ux - It is only when A/flux increases to about 10 6 we begin to see a slow 
down of field decay due to the presence of a relatively long period of forward creeping phase 
before the star eventually enters the co-moving phase. 

The magnetic evolution is more dramatic if we take 10 13 G as the initial field. Fig. [F]b 
shows that when AWx = 1, the star is basically in the co-moving phase all the time. When 
A/fl ux = 10 3 , a delay in field decay is observed due to presence of a long period of forward 
creeping phase at early times. Finally, when A/" nux = 10 6 , the star switches to an initial 
exponential, followed by an eventual power law delay mode, indicating that the star has 
locked into the forward creeping — > co-moving — > reverse creeping pattern. 

In summary, we find that in all reasonable parameter range, the value of the initial 
magnetic field will be reduced to 1/e of its original value in less than 1 Myr time when 
A/fl U x ~ 1 or when -B star (0) < 10 12 G. If A/fl ux ~ 10 6 , the characteristic field decay time is 
less than 20 Myr provided that B stSLI (0) < 3 x 10 12 G or 7 ~ 0. The only way we obtain a 
characteristic decay time longer than 20 Myr is by using A/flux as high as 10 6 , 5 star (0) as high 
as 10 13 G, and 7 ~ (see Fig. ||) — a set of physically unlikely combination of parameters. 



-24- 



4.6. Observational Implications Of Strange Star Field Decay 

The presence of a large number of systematic and random errors in the observed pulsar 
sample together with the fact that kinematic age sometimes does not truly reflect the real 
age of a pulsar greatly complicate the analysis of pulsar magnetic field decay. Nonetheless, 
current statistical analysis and computer simulation of galactic pulsar distribution suggest 
that the characteristic field decay time is at least 20 Myr (|Bhattacharya et al. 1992] ; |Wijers 



et al. 1993 ; see also the recent review by Bhattacharya and Srinivasan 1995). 



Various authors have studied theoretically the possibility of a magnetic field decay in 
the core of a neutron star flSauls 1989| ; |Srinivasan et al. 1990| ; |Chau et al. 1992] ; |Ding et aT 



. However, Pethick and Sahrling (1995) pointed out recently that the field decay time 
in the inner crust of a neutron star is at least 80 Myr due to the slow diffusion of the field 
through the inner crust. Consequently, core magnetic field decay only results in transporting 
the field from the stellar core to the inner crust; the total stellar field remains unchanged in 
> 80 Myr. Therefore, the conventional hypothesis that pulsars are neutron stars is consistent 
with the observed field decay time of the star. 



On the other hand, suggested by our numerical finding in §[4.5|, the magnetic field of 



a strange star decays with a characteristic times < 20 Myr. In fact, Eq. (pi ) tells us 
that the characteristic time for the buoyancy force dominated field decay rb UO y equals 4.2 x 
10~ 3 A/fl ux Myr. Thus, Tb U o y is longer than 20 Myr only when A/fl ux > 2 x 10 7 — a value which 
is attainable only when the quark superconductor is extremely type I. Therefore, our finding 
is inconsistent with the proposition that pulsars are strange stars. 



5. Conclusions 

In summary, I consider the rotation of a multi-superconducting species object. The 
rotation of such an object is made possible by the formation of vortices similar to that of su- 
perfluid helium. This finding implies the existence of vortex bundles in the core of a rotating 
superconducting strange star. Because it is energetically favorable for the vortex bundles to 
pin to magnetic flux tubes, the rotational dynamics and the magnetic field evolution of a 
strange star are coupled. 

Because the nuclear crust of a strange star cannot sustain a strong magnetic stress, the 
magnetic field evolution of the star is dictated by dynamics of the flux tubes in the core. The 
core magnetic field evolution due to inter-pinning of magnetic flux tubes and vortex bundles, 
and clumping of flux tubes is computed numerically. I find that in all reasonable parameter 
range, the characteristic decay time of the magnetic field is < 20 Myr. This finding does not 
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agree with the hypothesis that pulsars are superconducting strange stars because current 
pulsar data strongly suggest that pulsar magnetic fields do not decay in 20 Myr. 

A number of interesting questions remains. The phase diagram of interacting multiple 
component superconducting species in rotation remains unclear. It is also interesting to 
study the effects of alignment torque on magnetic and rotational history of a strange star, 
and the possible collapse of the nuclear matter crust due to magnetic stress. I plan to report 
them in my future works. 

I would like to thank Mai Ruderman for bringing up this problem to my attention 
and K. Y. Ding for providing us her field decay program for neutron stars for reference. 
I would also like to thank H.-K. Lo, Geoff Ravenhall and Frank Wilczek for their useful 
discussions. The hospitality of Aspen Center for Physics is acknowledged where the early 
part of this work was performed in the summer of 1995. This work is supported by DOE 
grant DE-FG02-90ER40542. 

A. Solution To The Strongly Interacting Superconducting Species Problem 

We prove the following claim by explicitly construct a solution to the problem. For 
simplicity, we shall drop all the star superscripts over the masses and charges. 

Claim: Suppose there is a finite number of superconducting species. Then Eq. ( pi]) has a 
solution with K ^ (and hence it has infinitely many solutions) if and only if there exist 
two vectors vi = (Mi, Qi) and v 2 = (M 2 , Q2) among the (rrij, qj) together with two rational 
numbers ai, a 2 such that (rrij,qj) = Pji^i + Pj2^2, with {1, /3ji, ^2} being linear dependent 
over the set of all rational numbers Q for all j. In addition, ctifyi +a 2 /3.,' 2 are rational numbers 
for all j, and a\Q 2 7^ a 2 Qi- 

Proof: (=>) By replacing vectors Vj by Vi/l (i = 1, 2) where i is the least common multiple 
of the denominators of a\(3j\ + a 2 /5j 2 , then (3ji will be replaced by Cf3ji. In addition, we 
can replace «j by Aa, where A is the least common multiple of the denominators of and 
«i/3ji+a 2 /3j 2 . Then it is easy to verify that A«j and ^(ai/3ji+a 2 /3j 2 ) are all integers. Thus, we 
can always assume that «i, a 2 and ai/3ji + a 2( 5j 2 to be integers instead of rational numbers. 

Suppose vi and v 2 are linear dependent on each other, then clearly the ratio of rrij to qj 
are the same for all j. So, by choosing Nj = for all j, K = 1, and & v = rnjc/qj, Eq. (|21~1 ) is 
satisfied. Now, we consider the more interesting case when vx and v 2 are linear independent 



-26- 



vectors. Consider the equation 

\ / A/f, —D, \ / Klh \ 

(Al) 



Pi) \M 2 -Q 2 J I $„//ic 



The linear independence of vi and v 2 implies that M\Q 2 — M 2 Q\ 7^ 0, and the solution of 
the above equation is given by 



(A2) 



( K/h \ = 1 / Q1P2 - Q2P1 \ 

{ § v /hc ) M 2 Q 1 - MiQ 2 \ M X P 2 - M 2 Px ) 

Clearly, by choosing P\ = a\ and P 2 = a 2 , we can check that the K given by Eq. ( |A2|) is 
non-zero. 

Now for each j, it is easy to check that by choosing Nj = ot\f3ji + a 2 (3j 2 (which is an 
integer, and this is possible because {1, /3ji, ^2} is linear dependent over Q) together with 
0), given in Eq. (|A2|) , then Eq. (|2lD is satisfied. (Clearly, if Nj, K, $ w is a solution 
of Eq. (0), then so is XNj, \K, \$ v for any non-zero integer A. Thus, Eq. d2lf) has either 
no solution, or infinitely many solutions.) □ 

(<=) On the contrary, suppose that for any linear independent vectors Vi and v 2 chosen 
among (rrij, qj), and any rational numbers ot\ and a 2 , we can find an j such that (rrij, qj) = 
PjiVi+f3j2V2 with either (1) {1, fyi, ^2} being linear independent over Q; or (2) aixfiji+aiiftji 
is irrational; or (3) OL1Q2 = ot-iQi- Now we analyze these three cases one by one: 



Case (1): If {1, fyx, (3j 2 } is linear independent over Q, then from Eq. (|A1| ), we know 
that Nj = P\f3ji + P2/3j2- In order that Pi, P 2 and Nj are all integers, the only possibility is 
Pi = P2 = Nj = 0. However, from Eq. (|A2|) , this implies K = and hence Eq. (|21|) has no 
solution for K ^ 0. 

Case (2): If ai/S^i + ct20j2 i s irrational, then similar to the argument in case (1), Nj = 
+ 012^2 is not an integer. Hence, this choice of ot\ and a 2 does not produce a solution 
for Eq. (P). 



Case (3): If (X1Q2 = ol 2 Q\, then K — from Eq. (|A2|) , which is impossible. □ 



The above claim implies that there is, in general, no solution to the strongly interacting 
limit when the number of different species of superconducting Cooper pairs is greater than 
or equal to three. On the other hand, solution of Eq. (plD always exists when there are only 
two species of Cooper pairs, as in the case of strange matter. 
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F/V A 




B 

Figl 

Fig. 1. — Schematic plot of F/V as a function of B for a fixed fl in a multi-superconducting 
non-interacting species sample. Dash lines represent the contributions from individual su- 
perconducting species, and the solid line is their combined free energy per unit volume. 
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Fig. 2. — The magnetic field decay of a strange star as a function of its initial B field. The 
initial field measured in G is varied from 10 13 (solid line), 3 x 10 12 (dash line), 10 12 (dotted 
line), 3 x 10 11 (dash dotted line), to 10 11 (short dash dotted line), (a) shows the field decay 
when A/fiux — 1; while (b) shows the field decay when A/fl ux = 10 6 . 
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Fig. 3. — Spin down of a strange pulsar as a function of its initial magnetic field. The solid, 
dotted, dash, dash dotted, and short dash dotted lines correspond to initial magnetic fields 
of 10 11 , 3 x 10 11 , 10 12 , 3 x 10 12 and 10 13 G respectively, (a) and (b) are for A/" flux = 1 and 
10 6 respectively. 
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Fig. 5. — Magnetic field evolution as a function of initial angular speed. The solid, dash and 
dotted lines correspond to initial angular speed of 200, 2000, and 12000 rads -1 respectively. 
A/fl U x is set to 10 6 . 
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Fig. 7. — Magnetic field evolution as a function of A/fl ux - The solid, dash, and dotted 
lines correspond to A/Vx = 1, 10 3 and 10 6 respectively, (a) shows the field decay when 
B(0) = 10 12 G, and (b) shows the field decay when B(0) = 10 13 G. Note that the solid curve 
overlaps with the dash one in (a). 
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